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ABSTRACT 

Young isolated radio-quiet neutron stars are still hot enough to be detectable at X-ray and 
optical wavelengths due to their thermal emission and can hence probe cooling curves. An 
identification of their birth sites can constrain their age. 

For that reason we try to identify the parent associations for four of the so-called Magnifi- 
cent Seven neutron stars for which proper motion and distance estimates are available. We 
are tracing back in time each neutron star and possible birth association centre to find close 
encounters. The associated time of the encounter expresses the kinematic age of the neutron 
star which can be compared to its characteristic spin-down age. Owing to observational un- 
certainties in the input data, we use Monte-Carlo simulations and evaluate the outcome of our 
calculations statistically. 

RX J 1856.5-3754 most probably originated from the Upper Scorpius association about 
0.3 Myr ago. RX J0720.4-3125 was either born in the young local association TW Hydrae 
about 0.4 Myr ago or in Trumpler 10 0.5 Myr in the past. Also RX J1605.3+3249 and RBS 
1223 seem to come from a nearby young association such as the Scorpius-Centraurus com- 
plex or the extended Corona- Australis association. For RBS 1223 also a birth in Scutum OB2 
is possible. 

We also give constraints on the observables as well as on the radial velocity of the neutron 
star. Given the birth association, its age and the flight time of the neutron star, we estimate the 
mass of the progenitor star. 

Some of the potential supernovae were located very nearby (< 100 pc) and thus should have 
contributed to the 10 Be and 60 Fe material found in the Earth's crust. 

In addition we reinvestigate the previously suggested neutron star/ runaway pair PSR 
B 1929+10/ C Ophiuchi and conclude that it is very likely that both objects were ejected during 
the same supernova event. 

Key words: stars: early-type - stars: kinematics - stars: neutron - pulsars: individual: PSR 
B1929+10, RX J1856.5-3754, RX J0720.4-3125, RX J1605.3+3249, RBS 1223. 



1 INTRODUCTION 

Neutron stars have very large proper motions which, with known 
distances, indicate high space velocities (e.g. Lyne & Lorimer 
1994; Lorimer et al. 1997; Hansen & Phinney 1997; Cordes & 
Chernoff 1998; Arzoumanian et al. 2002; Hobbs et al. 2005, there- 
after Ho05). Those high velocities that are usually larger than those 
of the progenitor stars are the result of asymmetric supernova ex- 
plosions assigning the new-born neutron star a kick velocity for 
that a number of mechanisms have been suggested (e.g. Burrows & 
Hayes 1996; Janka & Mueller 1996; Janka et al. 2005; Wang et al. 
2006; Kisslinger et al. 2009). 

Another fact is that about 46 % of O and 10 % of B stars also show 
high space velocities (see Stone 1991 for a discussion). Two sce- 
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narios (Hoogerwerf et al. 2001, thereafter H01) are accepted to pro- 
duce those so-called "runaway stars" (Blaauw 1961). The binary- 
supernova scenario is related to the formation of the high veloc- 
ity neutron stars: The runaway and neutron star are the products 
of a supernova within a binary system. The velocity of the former 
secondary is comparable to its original orbital velocity. The sec- 
ond scenario, which will not be further discussed within this paper, 
is the dynamical ejection due to gravitational interactions between 
massive stars in dense clusters. H01 proposed examples for each of 
the scenarios. In section 4 of this paper we will review and further 
investigate the binary supernova scenario for PSR B 1929+10 and 
the runaway 09.5V star ( Oph. Instead of restricting the unknown 
pulsar radial velocity to a certain interval, we will use a velocity 
distribution to cover the whole spectrum of possible radial veloci- 
ties. We then use more recent data for the pulsar to strengthen the 
hypothesis. 
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Beforehand, we introduce our sample of OB associations and clus- 
ters in section 2 and the procedure utilised in section 3. 
In section 5 we investigate the origin of four isolated radio-quiet 
X-ray emitting neutron stars. So far, seven such sources have been 
confirmed. They are not only bright X-ray sources, but are also de- 
tected in the optical due to a hot cooling surface. In two cases a par- 
allax is obtained. Brightness, parallax and temperature yield their 
radii. From spectra, one can in principle determine their mass and 
composition, which eventually may lead to constraints on the equa- 
tion of state. However, only seven such sources have been identified 
up to now for which they were named "The Magnificent Seven" 
(M7) (Treves et al. 2001; for a recent review see Kaplan 2008; two 
new candidates were recently published by Pires et al. 2009 and 
Rutledge et al. 2008). With known luminosity, one only needs to es- 
timate the age, in order to probe cooling curves. For young neutron 
stars, the characteristic spin-down age only gives an upper limit. 
We will therefore try to estimate their age from kinematics. 
Our aim in this paper is to identify the parent associations of four 
members of the M7 and find constraints on their observables and 
radial velocity. 

We give a summary of our results and draw our conclusions in sec- 
tion 6. 



2 THE SAMPLE OF ASSOCIATIONS AND CLUSTERS 

Given a distance of 1 kpc and typical neutron star velocities of 
100 to 500km/s (Arzoumanian et al. 2002; Ho05) and maximum 
ages of 5 Myr for neutron stars to be detectable in the optical (see 
cooling curves in Gusakov et al. 2005 and Popov et al. 2006), we 
restricted our search for birth associations and clusters of young 
nearby neutron stars to within 3 kpc. We chose a sample of OB 
associations and young clusters (we use the term "association" for 
both in the following) within 3 kpc from the sun with available 
kinematic data and distance. We collected those from Dambis et al. 
(2001) and H01 and associations to which stars from the galac- 
tic O-star catalogue from Mafz-Apellaniz et al. (2004) are associ- 
ated with. Furthermore, we added young local associations (YLA) 
from Fernandez et al. (2008) (thereafter F08) since they are possi- 
ble hosts of a few supernovae in the near past. We also included the 
Hercules-Lyrae association (Her-Lyr) and the Pleiades and massive 
star forming regions (Reipurth 2008a,b). We set the lower limit of 
the association age to ~ 2 Myr to account for the minimum life- 
time of a progenitor star that can produce a neutron star (progeni- 
tor mass smaller than « 30 M ; see e.g. Heger et al. 2003). 1 The 
list of all explored associations and their properties can be found 
in appendix A. Coordinates as well as heliocentric velocity com- 
ponents are given for a right-handed coordinate system with the x 
axis pointing towards the galactic centre and y is positive in the 
direction of galactic rotation. 



3 PROCEDURE 

To identify potential parent associations of neutron stars, we cal- 
culate the trajectories for both, the neutron star and the centre of 
the association, into the past. We account for solar motion adopting 

1 The minimum lifetime is actually as 6 Myr. We relax that condition to 
account for possible uncertainties in the association age as well as the fact 
that also more massive stars may produce neutron stars (Belczynski & Taam 
2008). 



a local standard of rest of (U V W) Q = (10.00 5.25 7.17) km/s 
(Dehnen & Binney 1998). To include the effect of a potential on the 
vertical motion we use the vertical acceleration from Perrot & Gre- 
nier (2003) (and references therein) and work with an Euler-Cauchy 
numerical method with a fixed time step of 10 4 years. Utilising this 
simple technique is fully sufficient for the treatment of some mil- 
lion years as done here and is consistent with results obtained by 
applying a fifth-order Runge-Kutta integration method for a more 
complicated potential as e.g. given in Harding et al. (2001) and ref- 
erences therein (cf. comparison of methods in Tetzlaff 2009). 
Simultaneously to the trajectories, we derived the separation be- 
tween the neutron star and the association centre for each time step 
and found the minimum of those as well as the associated time in 
the past. Since the input consists of observables with errors, we 
calculate a large set of simulations varying the starting parameters 
normally (proper motion and parallax 2 ; radial velocity see below) 
within their confidence levels. 

For neutron stars the radial velocity cannot easily be derived from 
spectra due to the large gravitational redshift and hence is un- 
known. To be able to investigate their three dimensional motion 
nonetheless, we varied the radial velocity within a probability dis- 
tribution derived from the one for three-dimensional pulsar ve- 
locities obtained by Ho05. They deducted a distribution with one 
Maxwellian component from transverse velocities of 233 pulsars 
assuming the radial component is of the same order. Arzoumanian 
et al. (2002) derived a distribution that includes two Maxwellian 
components suggesting that there might be two different popula- 
tions: those which gained their speed only due to a kick during 
the supernova and those which have an additional component from 
their former orbital velocity. However, as stated by Ho05 a one- 
component model better fits the observations. Infact, which radial 
velocity distribution is used is of minor importance since its usage 
shall just assure that the whole spectra of possible radial velocities 
is covered and the probability of occurence of excessively high val- 
ues is low which is the case for both distributions. Adopting the 
two-component model does not affect the final results significantly 
(Tetzlaff 2009). As the sign of the radial velocity cannot be derived 
from statistics it was chosen randomly, thus covering a total range 
of -1,500 to +1,500 km/s. 



4 A POTENTIAL BINARY SUPERNOVA IN UPPER 
SCORPIUS 

The runaway O-star £ Ophiuchi (£ Oph = HIP 81377) is an isolated 
main sequence star with a space velocity of ~ 20 km/s. Blaauw 
(1952) suggested its origin in the Scorpius OB2 association due to 
its proper motion vector which points away from that association. 
H01 investigated the origin of £ Oph in more detail and proposed 
that it gained its high velocity in a binary supernova in Upper 
Scorpius (US) about 1 Myr ago which is also supported by its 
large rotational velocity of > 340 km/s (Herrero et al. 1992; 
Penny 1996; Howarth & Smith 2001) and high helium abundance 
(Herrero et al. 1992; van Rensbergen et al. 1996; Villamariz & 
Herrero 2005). They also identified a neutron star which might 
have been the former primary: PSR B1929+10 (PSR J1932+1059). 

Initially, we repeat the experiment of H01 applying the 
same starting parameters for the pulsar (Table 1, /i* is the proper 

2 For some cases we obtain the parallax from other distance estimates. 



Identifying birth places of young isolated neutron stars 3 



motion in right ascension corrected for declination). Note that the 
pulsar radial velocity of 200 ± 50 km/s is not measured, but an 
assumption in HOI. The runaway proper motion and parallax were 
adopted from the Hipparcos Catalogue (Perryman et al. 1997) and 
its radial velocity from the Hipparcos Input Catalogue (Turon et al. 
1992). The results are in very good agreement with those of HOI. 
29,700 of 3 million Monte-Carlo runs yield a minimum separation 
between £ Oph and PSR B 1929+10 of less than 10 pc compared 
to 30,822 runs as published by H01. The smallest separation we 
found was 0.20 pc (0.35 pc found by H01). In 4,631 runs both 
objects were additionally not farther from the centre of US than 
10 pc. The latter number differs slightly from that of H01 (4,210 
runs) which is owing to somewhat different input parameters for 
the association. 

Similar calculations were also done by Bobylev (2008) with the 
same results using the epicycle approximation (Lindblad 1959; 
Wielen 1982) for tracing back the objects. Afterwards, he used 
more recent parameter values for the pulsar from Chatterjee et al. 
(2004) but increased the error intervals (factor of ten for the 
parallax, factor of 30 for the proper motion components), to have 
them of the same order as H01 before, and drew the conclusion 
that the scenario of a binary supernova in US involving the two 
objects is very likely. 

We reinvestigate the issue as well using a radial velocity dis- 
tribution (section 3) instead of a specified interval. In 3 million 
Monte-Carlo runs the smallest separation found was 4.0 pc. This is 
too large to support the hypothesis that ( Oph and PSR B1929+10 
once were at the same position at the same time in the past. How- 
ever, the errors of the parallax and the proper motion components 
of the pulsar (Table 1) are very small and thus might be underes- 
timated due to unknown systematical effects. For that reason we 
as well increase them by a factor of ten. The smallest difference 
in position then was 0.2 pc which is consistent with the binary su- 
pernova scenario. We select those minimum separations for which 
both, the pulsar and the runaway, were not more than 15 pc away 
from the centre of US. Of 3 million runs 5,367 fulfilled this re- 
quirement. The slope of the resulting histogram in Fig. 1(a) can 
be explained with a three-dimensional Gaussian distribution of the 
separation (see H01 for details), 
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respectively, either with y, = ± 5.6 pc or y, = 9.0 ± 4.0 pc. 
A denotes the three-dimensional separation \x pu i 3ar — x^o v h\ 
between the pulsar and the runaway, fi and a are the expectation 
value and the standard deviation, respectively, and tt = 3.1459 .... 
The agreement to the curves is remarkable. Fig. 2 shows the 
correlation between minimum separations and corresponding 
times. It can be seen that there is a small spread in time for larger 
separations which is owing to the deviations from a Gaussian 
distribution and the actual four dimensional problem since the time 
is involved as well. 

Adopting that PSR B1929+10 and £ Oph were ejected in the 
same supernova about 1 Myr ago in US, we find that the pulsar 




Figure 1. (a) Distribution of minimum separations d m i n of the 5, 367 runs 
for which both objects were not farther than 10 pc from the US centre 
with updated pulsar data. Drawn as well are theoretical curves for three- 
dimensional Gaussian distributions (equations 1 and 2) with /i = 9.0 pc 
and a = 4.0 pc (solid) and fi = and a = 5.6 pc (dashed), respectively. 
Note that these are not fitted to the data, (b) Distribution of corresponding 
flight times r in the past since the supernova. 




40 

(bin size: 2pc x 0.02Myr) 



Figure 2. r-d m i n contour plot to Fig. 1. 

would now have a radial velocity of ~ 250 km/s which is close to 
the value of 200 ± 50 km/s H01 adopted previously. The distance 
to the Sun of this supernova event 1 Myr ago was 154 ± 7 pc (note 
that this error quotes a 68% confidence interval but is not a la- 
error, though). 

Given an age of 5 Myr for US (de Geus et al. 1989; Preibisch 
et al. 2002) we can estimate the mass of the progenitor star of 
PSR B 1929+10 assuming that star formation happened contem- 
poraneously, i.e that the progenitor exploded (5 — 1) Myr after 
birth. Evolutionary models from Tinsley (1980) (thereafter T80), 
Maeder & Meynet (1989) (thereafter MM89) and Kodama (1997) 
(thereafter K97) (see also Romano et al. (2005) for a review) yield 
a progenitor mass of 44 Mq, 38 Mq and 32 Mq, respectively, 
for a star with a lifetime of 4 Myr assuming solar metallicity, 
corresponding to a spectral type of 06 to 07 on the main sequence 
(Schmidt-Kaler 1982). According to Sartori et al. (2003) US is 
somewhat older, namely 8 to 10 Myr. This association age range 
would imply a smaller progenitor mass of 18 (B0V) to 37 Mq 
(06V) (depending on the evolutionary model) which is in better 
agreement to the upper star mass limit of ~ 30 Mq (Heger et al. 
2003) for formation of a neutron star. However, owing to mass 
transfer during binary evolution also more massive stars may 



4 N.Tetzlaffetal. 



Table 1. Properties of PSR B 1929+10 and £ Oph (//* is the proper motion in right ascension corrected for declination). Note that the radial velocity of the 
pulsar is an assumption from HOI, not a measurement. 





«[°] 


8[°] 


7r [mas] 


/i* [mas/yr] 


US [mas/yr] 


v r [km/s] 


Ref. 


PSRB1929+10(H01) 


293.06 


10.99 


4±2 


99±12 


39±8 


200±50 


1 


PSR B 1929+ 10 (recent) 


293.06 


10.99 


2.76±0.14 


94.03±0.14 


43.37±0.29 




2 


COph 


249.29 


-10.57 


7.12±0.71 


13.07±0.85 


25.44±0.72 


-9.0±5.5 


3,4 



1 - H01 and references therein, 2 - Chatterjee et al. (2004), 3 - Hipparcos Catalogue, 4 - Hipparcos Input Catalogue 



produce neutron stars (Belczynski & Taam 2008). 
The 09.5V star £ Oph has a mass of 20 Mq (derived from 
isochrones by Schaller et al. (1992), Bertelli et al. (1994) and 
Claret (2004) for solar metallicity; table in Schmidt-Kaler 1982). 
If £ Oph was ejected in a supernova, the progenitor star should 
have an earlier spectral type, consistent with our results. With the 
additional constraint that the supernova produced a neutron star, 
i.e. the progenitor was later than «06 to 07, the progenitor star 
was 06/07 to 09. 

It may not be justified to enlarge the errors of the input param- 
eters of the pulsar. With the published small error bars we cannot 
confirm that £ Oph was at the same position in space at the same 
time as PSR B1929+10, but may have been ejected by a different 
supernova event in US. 



5 IDENTIFYING PARENT ASSOCIATIONS AND 
CLUSTERS FOR YOUNG ISOLATED NEUTRON 
STARS 

For four of the M7 distance estimates and proper motion values 
are available and thus their birth sites can in principle be found if 
assumptions on the radial velocity are made. Table 2 gives their 
properties. 

We trace back the objects for 5 Myr. As a first step we perform 
100,000 runs for each neutron star and all 140 associations in the 
sample. Those associations for which the smallest separation found 
was less than three times the association radius or less than 100 pc 
were chosen for a more detailed investigation with 2 million runs. 
We give the results of the latter in the following sections. 

5.1 RX J1856.5-3754 

Investigating probable birth places of RX J 1856.5-3754, nine 
of the 140 associations yield a smallest separation between the 
neutron star and the association centre found after another 2 
million runs which was consistent with the respective association 
radius 3 (see Table 3). However, only five of the other associations 
show a peak in the r-d min contour plot which also lies within 
the association boundaries or at least intersects. Table 4 lists the 
respective associations along with the results of the investigation 
including 2 million runs with varying input parameters. Columns 2 
and 3 of that table specify a region in the T-d min contour diagram 
with higher probability than its surroundings (see appendix B). 
Columns 4 to 8 give values for the present-day properties which 
RX J1856.5-3754 would have if it originated from the particular 

3 Considering the present radius of an association ensures that we do not 
miss potential supernovae that might have occured near the edge. The for- 
mer radius was surely smaller than the present one. 



Table 3. Associations for which the smallest separation between 
RX J1856. 5-3754 and the association centre found (min. d m i n ) after 2 mil- 
lion runs was within the radius of the association. Column 2 gives the small- 
est separation, column 3 the radius of the association (see Appendix A for 
references). 



Association 


min. d min [pc] 


RAssoc [PC] 


US 


0.4 


15 


Tuc-Hor 


46.1 


50 


P Pic-Cap 


26.3 


57 


Ext. R CrA 


4.4 


31 


AB Dor 


32.0 


45 


Her-Lyr 


10.4 


13 


Sgr OB5 


64.1 


111 


Sco OB4 


19.0 


33 


Pismis 24 


1.6 


2 


Tr27 


1.8 


6 



association. The space velocity w space which corresponds to the 
ejection speed of the neutron star in the supernova is given in 
column 7. Assuming that the ejection speed is of the same amount 
as the kick velocity, this value gives an upper limit of the latter. 4 
The last three columns indicate the position of the potential 
supernova at the corresponding time with the distance to the Sun 
and equatorial coordinates. For how those values are derived, 
please see appendix B. 

From Table 4 we conclude that the most likely parent asso- 
ciation for RX J1856.5-3754 is either Upper Scorpius (US), the 
extended Corona-Australis association (Ext. R CrA) or Scorpius 
OB4 (Sco OB4). For the two remaining - the 13 Pictoris group (J3 
Pic-Cap) and AB Doradus (AB Dor) - the radial needed velocities 
would be excessively high (> 1000 km/s, space velocities 
Vspace > 1000 km/s), thus making them unlikely (see Ho05) 
to host the birth place of RX J1856.5-3754. The probability of 
those two associations of being the birth associations is more than 
4, 000 times smaller than for the other three (for deviation of the 
probability, please see appendix B). 

Since the smallest separation found for Sco OB4 (19 pc) 
is relatively high compared to those of US (< 1 pc) and Ext. R 
CrA (< 5pc), and the modulus of the radial velocity needed of 
w 600 km/s (column 4 of Table 4) is rather high (space velocity 
w 665 km/s), Sco OB4 may be less likely the birth place of 
RX J1856.5-3754. For Ext. R CrA the potential place of the super- 
nova lies between 26 to 44 pc away from the association centre. 



4 The kick velocity strongly depends on whether the progenitor star was 
single or not and in the case of binarity also on the separation between the 
stars. However, supernova kick theory predicts kick velocity distributions 
which are consistent with observed neutron star velocities (e.g. Scheck et al. 
2006; Ng & Romani 2007; Kuranov et al. 2009). 
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Table 2. Parameters of four of the M7. Note that only for RX J1856. 5-3754 and RX J0720.4-3125 7r denotes trigonometric parallax measurements whereas for 
RX J1605.3+3249 and RBS 1223 distance limits are model estimates (from brightness and extinction compared to RX J1856.5-3754 and RX J0720.4-3125). 
References: 1 - Motch et al. (1999), 2 - Kaplan et al. (2002), 3 - Kaplan (2005), 4 - Schwope et al. (2005), 5 - Zane et al. (2006), 6 - Schwope et al. (2007) 
(their model 3), V. Hambaryan (priv. comm.), 7 - Posselt et al. (2007), 8 - Kaplan et al. (2007), 9 - Motch et al. (2007), 10 - Hohle et al. (2009). 







8[°] 


7r [mas] 


/i* [mas/yr] 


ti s [mas/yr] 


n H [10 20 cm" 2 ] 


Ref. 


RXJ1856.5-3754 


284.15 


-37.91 


5.6±0.6 


326.7±0.8 


-59.1±0.7 


0.74±0.10 


2, 3,7 


RXJ0720.4-3125 


110.10 


-31.43 


2.77±1.29 


-93.9±2.2 


52.8±2.3 


1.04±0.02 


8, 10 


RXJ1605.3+3249 


231.33 


32.82 


>2.4 


-43.7±1.7 


148.7±2.6 


1.1±0.4 


1,5,9 


RBS 1223 


197.20 


21.45 


1.4... 13 


-207±20 


84±20 


1.8±0.2 


4, 6,9 



Table 4. Potential parent associations of RX J1856. 5-3754. 

Columns 2 and 3 mark the boundaries of a 68% area in the T-d m i n contour plot for which the current neutron star parameters (columns 4 to 7, radial velocity 
v r , proper motion ^t* and fig and parallax 7r) were obtained and columns 8 to 10 indicate the distance to the Sun and equatorial coordinates (J2000.0) of 
the potential supernova. Column 7 gives the space velocity (ejection speed) v spaC e derived from proper motion and radial velocity. For the deduction of the 
values given in columns 4 to 11, please see appendix B. 



Association d min r v r /i* fig v space it d,Q a 8 



[pc] [Myr] [km/s] [mas/yr] [mas/yr] [km/s] [mas] [pc] [°] [°] 



US 


5. 


.13 


0.28. . 


.0.34 


193t« 


326.7±0.8 


-59.1±0.7 


349l 4 3 ° 


5 41+0.33 


140 .. . 163 


243.211°;^ 


-23.66±«;H 


13 Pic-Cap 


35. . 


.53 


0.10. . 


.0.15 


1190±1£ 


326.7±0.8 


-59.1±0.7 


1222^^3 


r: 71 +0.49 
°- (1 -0.24 


40t* 


221... 233 


-18. . . - 8 


Ext. R CrA 


26. 


.54 


0.10. . 


.0.15 


406^ 21 


326.7±0.8 


-59.1±0.7 


492l^ 3 


5 68+ ' 35 


127^ 


256.76l 2 ;^ 


-33.721°;^ 


AB Dor 


44. 


.60 


0.12. . 


.0.16 


1146111 


326.7±0.8 


-59.1±0.7 


1177 -m 


c oo+0. 30 
5.88_ 2? 


32... 46 


210.19l|;^ 


0.90l 2 ;^ 


Sco OB4 


25. 


.33 


1.26. . 


.1.51 


-604±i£ 


326.7±0.8 


-59.1±0.7 


664l 2 ^ 


r 79 +0.32 
''■ ' -0.30 


1092± 2 


259.84t°;«| 


-31.66±»;lt 



The radius of Ext. R CrA is 31 pc, thus the supernova should have 
occurred near the edge of the association. In comparison to that, 
a supernova in US which produced RX Jl 856.5-3754 would have 
taken place at 5.5 ± 1.2 pc from the centre of US (Fig. 3) which 
has a radius of 15 pc. As supernovae are supposed to occur in 
denser regions of an association, the latter scenario thus seems 
most likely. Furthermore, the probability of US being the parent 
association of RX J1856. 5-3754 is 10 times larger than that for 
Ext. R CrA and 8 times larger than for Sco OB4. 
The origin of RX J 1856.5-3754 has also been previously suggested 
by Walter & Lattimer (2002) to lie within the US association 
~ 0.5 Myr ago as the projected trajectory of RX J1856. 5-3754 
crosses the projection of US. Here, we find a kinematic age of 
w 0.3 Myr. 

In this case (Fig. 4) the neutron star would now have a radial 
velocity of £s 190 km/s. The space velocity of ^ 340 km/s in 
this case is close to the peak of the velocity distribution for pulsars 
from Ho05. 

Walter et al. (2000) suggested that ( Oph was the former 
secondary of RX J1856. 5-3754 that exploded a few Myr ago in US. 
Therefore, we re-investigate the separation between RX J1856.5- 
3754 and the runaway star as well as their separation to the centre 
of US. Out of 3 million runs the smallest separation found between 
the neutron star and £ Oph was 22 pc which is far too large to 
support this former binary hypothesis. 

Given the age of US of 5 Myr (de Geus et al. 1989) we 
estimate the mass of the progenitor star (lifetime 4.7 Myr) to 
be, depending on the model, 45 M Q (T80), 41 M (MM89) and 
35 M Q (K97) (spectral type ^06.5, Schmidt-Kaler 1982). Using 
the more recent age determination for the Sco-Cen associations by 
Sartori et al. (2003) of 8 to 10 Myr for US, we estimate the mass 
of the progenitor star of RX 11856.5-3754 to be between 17 and 



35 M (B0 to 07 on the main sequence). Currently, the earliest 
spectral type in US is B0. Hence, our assumption of contemporary 
star formation and our result on the progenitor spectral type are 
consistent with the progenitor star of RX 11856.5-3754 being 
earlier than the earliest present member star. 

Since Walter & Lattimer (2002) published a significant differ- 
ent parallax for RX J1856. 5-3754 of n — 8.5 ± 0.9 mas, we repeat 
our simulations adopting this value. The overall picture of poten- 
tial parent associations remains (cf. Table 3 5 ). Still, US is the most 
probable association to have hosted the supernova which formed 
RX J 1856.5-3754. While the position of this supernova using Wal- 
ter's parallax does not significantly change, the neutron star's cur- 
rent parallax would be ty = 7.44l(,'4g mas and the radial veloc- 
ity needed v r = —25^24 km/s. Furthermore, the kinematic age 
would be slightly larger, r w 0.5 Myr (confirming the result of 
Walter & Lattimer 2002), implying a progenitor mass of 37 to 
45 Mq (~06.5 on the main sequence) adopting an age of US of 
5 Myr or 17 to 29 M (B0 to 07 on the main sequence) for an as- 
sociation age of 10 Myr (T80, MM89, K97). 
We again investigated the RX 11856.5-3754/ ( Oph binary scenario 
using the larger parallax as an input. Out of 2 million runs we found 
a smallest separation of 20 pc ruling out this hypothesis. 

5.2 RX J0720.4-3125 

For RX 10720.4-3125, 12 associations were identified for which 
the smallest separation between the neutron star and the associa- 
tion centre found in the past after 2 million runs lied within the 
respective association boundaries (see Table 5). Five of those show 
a peak in the T-d min contour plot at values of d m in that at least 

5 We find additional small separations for Bochum 13 and NGC 6383; 
however, the fraction of such runs is very small. 
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Figure 3. Distribution of minimum separations d min between RX J1856.5- 
3754 and the centre of US within the defined time range r (see Table 4, 
column 3). Shown as well is a theoretical curve for a three-dimensional 
Gaussian distribution (equation 1) with fj, = 5.5 pc and a = 1.2 pc. Note 
that the curve is not a fit to the data but shall just give an explanation of the 
slope of the histogram. 

The disagreement on the right side of the diagram is due to the actual four 
dimensional problem since the time is involved. Furthermore, natural dis- 
crepancies to a Gaussian distribution occur since parallaxes are taken to 
obtain positions and a Maxwellian distribution is utilised for the velocities. 
Unlike in the previous section we cannot restrict to a certain number of runs 
owing to a third component (in the case of PSR B 1929/ C Oph the distance 
to US) which gives constraints on the time. 



+1E°- 




: 0" 



.,■..301 



P! 330° 



Figure 4. Past trajectories for RX J1856.5-3754 and US projected on a 
Galactic coordinate system (for a particular set of input parameters consis- 
tent with Table 4). Present positions are marked with a star for the neutron 
star and a diamond for the association. The large circle reflects an associa- 
tion radius of 15 pc. 

intersect the association. For the Chamaeleon-T association (Cha 
T) only a very small fraction of runs yield small distances, thus we 
exclude it as a probable parent association for RX J0720.4-3125. 
Analogues to Table 4 for RX J1856.5-3754, they are given in 
Table 6 along with the separation and corresponding time ranges 
in the past defined by the peak in the T-d m i n contour diagram. 
Also given are the neutron star parameters needed as well as the 
position of the potential supernova. 



Table 5. Associations for which the smallest separation between RX 
J0720.4-3125 and the association centre found (min. d min ) after 2 million 
runs was within the radius of the association, as Table 3. 



Association 


min. d min [pc] 


RAssoc [PC] 


US 


11.8 


15 


UCL 


27.5 


33 


LCC 


18.2 


23 


TWA 


0.1 


33 


Tuc-Hor 


37.6 


50 


/3 Pic-Cap 


18.2 


57 


HD 141569 


4.7 


16 


Pup OBI 


6.0 


65 


PupOB3 


0.5 


15 


NGC 2546 


6.6 


10 


Tr 10 


17.3 


23 


ChaT 


0.8 


21 



It has been previously suggested that RX J0720.4-3125 orig- 
inated from the Trumpler 10 (Tr 10) association by Motch et al. 
(2003) who considered the general direction of the neutron star's 
motion and Kaplan et al. (2007) who investigated the probability of 
close approaches of the neutron star to any of the OB associations 
given in de Zeeuw et al. (1999) giving the proper motion randomly 
directions and sign as well as varying the radial velocity in a 
certain range and adopting specified distances. They concluded 
that it is very likely that Tr 10 is the birth place. Since early 
B-type stars are present in Tr 10, it is plausible that the association 
already experienced a supernova. Varying the parallax in the range 
n ± and the radial velocity in the range v r ± 0.935«t (0.935 
corresponds to la in v r for random orientation, vt is the transverse 
velocity), Kaplan et al. (2007) found a separation between the 
neutron star and the centre of Tr 10 of 17 pc 0.7 Myr ago for a 
radial velocity of —20 to +50km/s. Infact, we find the smallest 
separation between RX J0720.4-3125 and the centre of Tr 10 to be 
17 pc ^ 0.6 Myr in the past. We find radial velocities in the range 
of -50 to 500 km/s for which RX J0720.4-3 125 would have come 
as close as 30 pc to the centre of Tr 10. In the T-d mi „ contour plot 
a maximum lies at d min between 23 and 40 pc and corresponding 
flight times r of 0.44 and 0.65 Myr. The slope of the histogram 
of minimum separations d min within this time range (Fig. 5(c)) 
is in good agreement with a three-dimensional Gaussian implying 
d m in = 33.0 ± 2.8 pc. Given a radius of Tr 10 of 23 pc, this 
means that the potential supernova would have occurred outside 
the association or, if we relax the boundary condition, at least at 
its outer edge making the scenario less likely (the probability is 
two times smaller than for TWA). If we nevertheless adopt Tr 10 
as the birth place of RX J0720.4-3125 (Fig. 6) and consider the 
association age (15 to 35 Myr), we estimate a progenitor mass of 
w 7 to 13 M Q (models from T80, MM89 and K97), corresponding 
to a spectral type between B4 and Bl on the main sequence 
(Schmidt-Kaler 1982). Motch et al. (2003) and Kaplan et al. 
(2007) also suggest Vela OB2 (Vel OB2) as a possible birth place 
of RX J0720.4-3125 for which we find a minimum separation of 
81.6 pc (cf. 70 pc stated by Kaplan et al. 2007) which is outside 
the association with a radius of 35 pc. 

We exclude Pup OB3 as birth association owing to a very small 
parallax needed to reach the association implying a current 
distance to the Sun of > 1 kpc. For that reason, the probability of 
Pup OB3 being the birth site of RX J0720.4-3125 is more than five 
times smaller than for TWA (and more than two times for most of 
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Table 6. Potential parent associations of RX J0720.4-3125, columns as in Table 4. 



Association 



[pc] 



[Myr] 



[km/s] [mas/yr] [mas/yr] [km/s] 



[mas] 



do 
[pc] 



a 
[°] 



S 

[°] 



TWA 


0. 


.12 


0.34. 


.0.46 


5021H 1 


-93.9±2.2 


52.8±2.3 




4 06+ ' 47 

4.UO_q 4 g 


48.. 


Tuc-Hor 


44. 


.80 


0.26. 


.0.52 


476± l 7 f 


-93.8±2.2 


52.8±2.3 


4921 


4 in+0.55 


35.. 


13 Pic-Cap 


35. 


.84 


0.31. 


.0.61 


468t^ 6 


-93.9±2.2 


52.8±2.3 


486l|° 7 


3.94lg;» 


38.. 


HD 141569 


12. 


.22 


0.52. 


.0.75 


469t^ 8 


-93.8±2.2 


53.0±2.3 




3.80t° j6 


96... 


Pup OB3 


14. 


.42 


0.47. 


.0.61 




-94.0±2.1 


52.9±2.3 


71 5+ 143 


0.84±°;°* 


1440 


Tr 10 


23. 


.40 


0.44. 


.0.65 


2go +143 


-93.8±2.2 


53.1±2.2 




i qo+0.24 


335.. 



.67 
.70 
.78 
116 

fl2 

-23 

.375 



187.10± 3 2 . 44 



190.19 



+4.00 
4.97 

209.87± 3 ; 2 3 l 
240 .. . 250 

1 24 1 1 + 1 - 27 
135 .. . 140 



-37.98+J;3g 
-37.29t 1 1 3 ° 
-28.52l 2 2 ;? 4 



-37.16l . 43 5 
-4D 1 q+ - 58 



the other associations listed in Table 6). 

The four other associations for which we find a maximum 
in the T-d m i n contour plot which is located well within the 
association boundaries are young local associations - TW Hydrae 
(TWA), Tucana/Horologium (Tuc-Hor), the /3 Pictoris moving 
group (J3 Pic-Cap) and HD 141569. For two of them, TWA 
and the HD 141569 group, the smallest separation between RX 
J0720.4-3125 and the centre of the particular association is smaller 
than 5 pc. For that reason it seems plausible that one of those two 
may be the birth place of RX J0720.4-3125. However, since the 
radii of Tuc-Hor (50 pc) and (5 Pic-Cap (56 pc) are rather large, 
they could also have hosted the supernova. Compared to TWA 
the probabilities that one of those two associations hosted the 
supernova which formed RX J0720.4-3125 is two to three times 
smaller. 

Figs. 5(a) and 5(b) show the distributions of minimum sep- 
arations d min to the centres of TWA and HD 141569. From the 
theoretical curves we can predict that a potential supernova which 
formed RX J0720.4-3125 in HD 141569 would have been situated 
15 ± 2 pc from the centre, i.e. at the edge of the association with a 
radius of about 15 pc. Moreover, since HD 141569 only consists of 
five members (F08), it is rather unlikely to have hosted a supernova 
far from its core. The probability of such a scenario is seven times 
smaller than for TWA. 

In the case of TWA we can predict the supernova very close to the 
association centre at a separation of ± 2 pc. Then, RX J0720.4- 
3 125 would have been formed ^0.4 Myr ago and its current radial 
velocity would be w 500 km/s (space velocity w 520 km/s). The 
current distance to the Sun which the neutron star would have in 
this case is w 250 pc which is in very good agreement with 235 to 
265 pc as derived by Posselt et al. (2007) (thereafter P07) from the 
spectrum and uh (note that they obtained a slightly higher value 
of 1.2 • 10 20 cm~ 2 compared to the most recent value from Hohle 
et al. 2009; Table 2) as well as with the parallax of 2.77± 1.29 mas 
(Kaplan et al. 2007). 

The potential supernova would have occurred between 48 and 
67 pc from the Sun. For such a small distance we expect to 
find supernova-produced radionuclides in the terrestrial crust (El- 
lis et al. 1996) such as 10 Be and 60 Fe, with half-lives of 1.5 Myr 
(Acton et al. 1993) and 2.6 Myr (Rugel et al. 2009), respectively. 
A small but insignificant signal of 60 Fe was found at a time below 
1 Myr (Knie et al. 2004; Fitoussi et al. 2008). For 10 Be anomalous 
abundances in Antarctic ice cores corresponding to times of « 0.35 
and « 0.6 Myr in the past were identified (Raisbeck et al. 1987). 
The first signal has also been reported for deep-sea sediments in the 
Gulf of California (McHargue et al. 1995) and in the Mediterranean 




(a) TWA (b) HD 141569 




15 20 25 30 35 40 

(c) Tr 10 



Figure 5. (a) Distribution of minimum separations d min between RX 
J0720.4-3125 and the centres of TWA (a), HD 141569 (b) and Tr 10 (c), re- 
spectively, within the defined time range r (see Table 6, column 3). Shown 
as well are theoretical curves for three-dimensional Gaussian distributions 
(equation 1) with fi = Opc and a = 2pc for TWA, fi = 15 pc and 
a = 2 pc for HD 141569 and jx =33pcand<r = 2.8 pc for Tr 10. Note 
that the curves are not fitted to the data (see also Fig. 3). 

(Cini Castagnoli et al. 1995). More recently, Cole et al. (2006) anal- 
ysed ice cores from the Qinghai-Tibetan plateau in China but could 
not confirm previous findings. However, it seems possible that a re- 
cent supernova in TWA which produced RX J0720.4-3125 (Fig. 6) 
contributed to one of the 10 Be signals. 

Moreover, F08 state that it is likely that one or more of the YLA 
hosted a supernova in the near past. Despite the small number of 
TWA members, we can expect one 10 M Q -star to have been born 
in TWA due to its present mass function 6 (Fig. 7). Given the asso- 
ciation age of 3 to 20 Myr (e.g. Barrado Y Navascues 2006), we 

6 We estimated the masses for the TWA members listed in F08, taking also 
into account that HD 98800 (Soderblom et al. 1998) and TWA-5 (Neuhauser 
et al. 2000; Torres et al. 2003) are both quadruple, using the D'Antona & 
Mazzitelli (1997) model, and 2.9 M for the AO-type main sequence star 
HR 4796 A. 
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Figure 6. Past trajectories for RX J0720.4-3125 and Tr 10 and TWA, re- 
spectively, projected on a Galactic coordinate system (for particular sets 
of input parameters consistent with Table 6). Present positions are marked 
with a star for the neutron star and a diamond for Tr 10 and an open circle 
for TWA. Large circles reflect association extensions (radii of 23 pc for Tr 
10 and 33 pc for TWA). 




mass [M syn ] 



Figure 7. Present TWA mass function (black circles, see footnote 5) in com- 
parison with the initial mass function £ (IMF, solid line) as given by Kroupa 
& Weidner (2005), their equation 2. Dashed and dotted lines represent the 
la and 1.5a IMF boundaries. Within 1.5a we expect up to one star with 
10 Mq ; the probability for higher mass stars is low. 

estimate a progenitor mass of 10 (B2V) up to 100 M (04V) (10 
to > 50 M for T80, 11 to 100 M Q for MM89 and 10 to 100 M 
for K97); spectral types from Schmidt-Kaler (1982)). 
If RX J0720.4-3125 was born in TWA, the supernova would have 
contributed only a small amount to 10 Be and 60 Fe found in the 
Earth's crust due to the relatively low mass progenitor (w 10 M©). 
This would be consistent with the findings. Furthermore, we pro- 
pose a larger age (« 20 Myr) for TWA. 

5.3 RX J1605.3+3249 

Since there is currently only an upper limit of the distance 
« 410 pc, Motch et al. 2007) available for RX J1605.3+3249, 
thus a lower limit for the parallax (n > 2.4 mas, see Table 2) we 
had to adopt a value for our investigations. The hydrogen column 
density of n H = 1.1 ± 0.4 • 10 20 cm~ 2 (Motch et al. 1999) 
predicts a relatively small distance of RX J1605. 3+3249 to the 
Sun. For that reason we chose a parallax value of n = 5 ± 3 mas 
corresponding to a distance of 200^5° pc. Infact, this is a 



Table 7. Associations for which the smallest separation between RX 
J1605. 3+3249 and the association centre found (min. d min ) after 2 mil- 
lion runs was within the radius of the association, as Table 3. 



Association 


min. dmi n [pc] 


RAssoc [PC] 


US 


4.6 


15 


UCL 


13.0 


33 


Tuc-Hor 


5.2 


50 


P Pic-Cap 


6.3 


57 


eCha 


21.4 


28 


HD 141569 


9.1 


16 


Ext. R CrA 


25.3 


31 


AB Dor 


12.7 


43 


NGC 6322 


2.2 


11 


Bochum 13 


2.1 


2 


Sco OB4 


0.9 


33 



source of large uncertainty. After the 2 million run investigation 
we identify 11 associations for which the smallest separation 
dmin between the association centre and RX J1605. 3+3249 
was consistent with the association radius (Table 7). Five of 
them show a maximum in the T-d m in contour plot also being 
consistent with the respective association boundaries (Table 8, top). 

Due to the small parallax needed for RX J1605.3+3249 to 
reach the Scorpius OB4 (Sco OB4) association (n — 1.1 mas, col- 
umn 7 of Table 8) we may exclude this association since this is even 
smaller than the lower limit of ty > 2.4 mas (Table 2). 
The remaining four associations once again belong to the YLA 
(F08). We can explain the slopes of the histograms of minimum 
separations d m in with three-dimensional Gaussians (equation 1) 
that predict the place of the potential supernova within the bound- 
aries of the associations (see Fig. 8). For the YLA, the simulations 
suggest a very small current distance of the neutron star to the Sun 
between 100 and 250 pc. From rin P07 derived a distance between 
325 and 390 pc (they obtained n H = 2.0 • 10 20 cmT 2 without 
giving errors; however, for an nn of « 1 • 10 20 cm~ 2 (Table 2), 
the distance should at least be w 230 pc, cf. P07 their Table 2). 
For that reason we find Tuc-Hor, j3 Pic-Cap and AB Dor to be less 
likely parent associations of RX J1605. 3+3249 since the neutron 
star would then need to have a current distance to the Sun smaller 
than 160 pc. 7 

In Table 8 (bottom) we also give the results for two associations of 
the Scorpius-Centaurus (Sco-Cen) complex - US and Upper Cen- 
taurus Lupus (UCL) - for which slightly larger distances are re- 
quired. Constraining the age to 1 Myr as a maximum and using 
current distances of 100 to 500 pc as well as radial velocities in 
the range of ±700 km/s Motch et al. (2005) proposed the origin of 
RX J1605.3+3249 in US. We consider the Sco-Cen complex also 
to be a likely birth site of RX J1605.3+3249. The respective d mm 
histograms are shown in Figs. 8(e) and 8(f). 

Hence, we conclude that RX J1605. 3+3249 either originated from 
the extended Corona-Australis association « 0.5 Myr ago with a 
current distance of « 230 pc and radial velocity of £s 320 km/s 
(space velocity « 360 km/s) or was born in the Sco-Cen complex 
28.5 ± 2.0 pc from the centre of US about 0.5 Myr in the past 
(Fig. 9). The first scenario is possible due to the existence of B- 
type stars (highest content within the YLA, F08), thus late O- or 

7 Estimating probabilities, even rough values, is unreliable in the case of 
RX J1605. 3+3249 because of the uncertain distance. Hence, an individual 
discussion is appropriate. 



Identifying birth places of young isolated neutron stars 9 



Table 8. Potential parent associations of RX J1605. 3+3349, columns as in Table 4. 



Association d min r v r /i* fig v space n cLq a S 

[pc] [Myr] [km/s] [mas/yr] [mas/yr] [km/s] [mas] [pc] [°] [°] 



Tuc-Hor 


7. 


.15 


0.24. 


.0.39 


452l^ B 


-43.7±1.7 


148.7±2.6 


462±^ 8 


7.54: 


hi. 53 
-1.33 


3911 


o ,,0 + 10.60 
32.4S_ 34 55 


-75.69: 


h4.70 
-1.55 


13 Pic-Cap 


8. 


.10 


0.16. 


.0.24 




-43.7±1.7 


148.6±2.6 


543^° 4 4 


8.86: 


-1.60 
-1.09 


18±? 


292.491^-51 


-75.21: 


h8.25 
-1.41 


Ext. R CrA 


30. 


.42 


0.43. 


.0.64 


32118" 


-43.7±1.7 


148.8±2.6 


362l^ 7 


4.41: 


1-0.90 
-0.85 


70... 115 


261.99lJ;^ 


-39.84: 


h2.10 
-2.00 


AB Dor 


21. 


.43 


0.17. 


.0.32 


452±|°° 


-43.7±1.7 


148.7±2.6 


462^11° 


7.56: 


hi. 22 
-0.78 


21+ 7 
zr_ 6 


264.30H™ 


-44.11: 


h6.98 
-8.82 


Sco OB4 


12. 


.74 


1.28. 


.1.68 




-43.7±1.7 


148.7±2.6 


7431^ 


i.09: 


h0.08 
-0.08 


1042... 1140 


259.571°;^ 


-35.0. . . - 


30.5 



US 26. ..36 0.46. ..0.63 343l|g -43.7±1.6 148.8±2.6 4O6+I5 3.38+^; 129 . . . 171 255.09+o;gg -18.44 

UCL 50. ..63 0.48. ..0.73 318l 6 3 6 -43.7±1.7 148.7±2.6 375l!,2 2 3.70l ;g9 96 . . . 164 260.48l l ;|g -35.68 




(a) Tuc-Hor 



(b) (3 Pic-Cap 




(c) Ext. R CrA 




(d) AB Dor 




(f) UCL 



Figure 8. (a) Distribution of minimum separations d min between RX 
J 1605. 3+3249 and the centres of Tuc-Hor (a), /3 Pic-Cap (b), Ext. R CrA 
(c), AB Dor (d), US (e) and UCL (f), respectively, within the defined time 
range r (see Table 8, column 3). Shown as well are theoretical curves for 
three-dimensional Gaussian distributions (equation 1) with fi = 9.5 pc 
and a = 1.0 pc for Tuc-Hor, ft = 8.5 pc and a = 0.4 pc for /3 Pic- 
Cap, = 33.5 pc and a = 1.5 pc for Ext. R CrA, = 29.0 pc and 
a- = 4.5 pc for AB Dor, fi = 28.5 pc and a = 2.0 pc for US and 
H = 53.5 pc and a = 2.5 pc for UCL. Note that the curves are not fitted 
to the data (see also Fig. 3). 




Figure 9. Past trajectories for RX J1605.3+3249 and US and Ext. R CrA, 
respectively, projected on a Galactic coordinate system (for particular sets 
of input parameters consistent with Table 8). Present positions are marked 
with a star for the neutron star and a diamond for US and an open circle for 
Ext. R CrA. Large circles reflect association extensions (radii of 15 pc for 
US and 31 pc for Ext. R CrA). 



early B-type stars may have been born there. For the latter scenario 
the current distance to the Sun would be w 300 pc and the neutron 
star's radial velocity w 340 km/s (space velocity w 400 km/s). 
Both cases would be in good agreement with a space velocity of 
w 400 km/s which reflects the mean velocity of the velocity dis- 
tribution for pulsars from Ho05. 

Accepting one of the two birth sites and taking the ages of the as- 
sociations of 5 to 10 Myr (US) and 10 to 15 Myr (Ext. R CrA), we 
estimate the mass of the progenitor star to be 17 (B0V) to 45 Mq 
(05.5V) (29 to 45 M for T80, 18 to 42 M for MM89, 17 to 
36 M for K97) for a birth in the Sco-Cen complex or 12 (B1V) to 
30 M (07V) (13 to 30 M for T80, 13 to 18 M for MM89, 12 
to 17 M Q for K97) for a birth in Ext. R CrA (spectral types from 
Schmidt-Kaler 1982). 

Such a nearby supernova (see Table 8, column 8) could also have 
contributed to the 10 Be and 60 Fe material found in the Earth's crust 
(e.g. Raisbeck et al. 1987; Knie et al. 2004). 



5.4 RBS 1223 (RX J1308.8+2127) 

For RBS 1223 both distance limits are available (76 . . . 700 pc, 
Schwope et al. 2005; Motch et al. 2007, hence ■k = 1.4 . . . 13 mas). 
Since the span between those limits is rather large, we adopted a 
parallax of tt — 7 ± 5 mas which is, as for RX J1605. 3+3249, 
a source of large uncertainty. After performing 2 million Monte- 
Carlo runs, we found 22 associations for which the smallest 
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Table 9. Associations for which the smallest separation between RBS 1223 
and the association centre found (min. d m ; n ) after 2 million runs was 
within the radius of the association, as Table 3. 



Association 


min. d mi „ [pc] 


R Assoc 


US 


1.1 


15 


Tnp-l-Tnr 

1 U-W .1 1U1 


10.3 


50 


R Pir-Can 

L V„ clIJ 


0.3 


57 




0.1 


16 


Fxt R OA 


0.3 


31 


AB Dor 


7.9 


43 


Hcr~Lyr 


9.7 


13 






in 


NGC 6530 


1.3 


6 


SgrOBl 


2.2 


104 


Sgr OB7 


3.0 


29 


Sgr OB4 


5.0 


27 


Sgr OB6 


1.8 


17 


M 17 


3.5 


10 


SerOBl 


1.1 


36 


NGC 6611 


0.8 


7 


Set OB3 


2.2 


17 


Ser OB2 


3.7 


31 


NGC 6604 


4.5 


8 


Set OB2 


6.6 


31 


NGC 6383 


3.0 


5 


M6 


2.5 


4 




d min [PC] d min [pc] 



(a) Ext. R CrA (b) US 

Figure 10. Distribution of minimum separations d min between RBS 1223 
and the centre of the extended Corona-Australis association (a) and US (b) 
within the defined time range r (see Table 10, column 3). Shown as well 
is a theoretical curve for a three-dimensional Gaussian distributions (equa- 
tions 1 and 2) with [m = Opc and a = 13.5 pc for Ext. R CrA and 
^ = 51 pc and a = 12 pc for US. Note that the curves are not fitted to 
the data but shall just give an explanation of the slopes of the histograms. 




- : +3Qf- 



separation between RBS 1223 and the association centre found 
was within the association radius (Table 9). Five of these 23 show 
well defined areas of higher probability in their T-d min contour 
plots (Table 10), whereas for more distant (> 1 kpc) associations 
as well as Serpens OBI (Ser OBI) and M6 (below Her-Lyr in 
Table 9) only very few runs yield small separations d min between 
the neutron star and the association centre. The latter behaviour 
does not change even if we constrain input parallaxes to a range of 
3 ± 2 mas. Furthermore, an origin in most of those associations 
implies a parallax smaller than the lower limit of 1.4 mas. For 
those reasons we may exclude them as probable birth associations 
for RBS 1223. 

However, calculating past trajectories of RBS 1223, Motch 
et al. (2007) proposed that a small fraction of orbits passes though 
the Scutum OB2 (Set OB2) association. More recently, Motch et al. 
(2009) reinvestigated this association considering that it consists 
of two groups, Set OB2A at 510 pc and Set OB2B at 1170 pc 
(Reichen et al. 1990). We adopt these two distances and the proper 
motion and radial velocity for Set OB2 from Dambis et al. (2001) 
and calculate again 2 million orbits for both, the neutron star and 
each subgroup, as well as their separations d m i n at each time step 
r. Due to the large uncertainty in parallax and radial velocity of 
RBS 1223, we do not find a maximum in their T-d m i„ plots, but 
can derive flight time ranges from the r histograms for which 
smaller separations {d min ^ 50 pc) occur. Considering the input 
parameters of the runs in the derived time ranges, we obtain the 
neutron star properties and potential supernova positions as given 
in Table 11. 

All associations for which we found a peak in their T-d min 
contour plot that was consistent with the association boundaries 
(Table 10, top) belong to the YLA (F08) and imply a rather small 
present distance of RBS 1223 to the Sun. Regarding the hydrogen 
column density of rin = 1.8 ± 0.2 • 10 20 cm~ 2 (Schwope et al. 




Figure 11. Past trajectories for RBS 1223 and US, Ext. R CrA and Set OB2, 
respectively, projected on a Galactic coordinate system (for a particular set 
of input parameters consistent with Table 10). Present positions are marked 
with a star for the neutron star, a diamond for US, an open circle for Ext. R 
CrA and a square for Set OB2. Large circles reflect association extensions 
(radii of 15 pc for US, 31 pc for Ext. R CrA) and 31 pc for Set OB2. 



2007) which is comparable to other M7 members and even a bit 
higher (see Table 2), the distance of RBS 1223 to the Sun is sup- 
posed to be at least 180 pc (the distance to the Sun of the probably 
closest member RX J1856. 5-3754 is « 180 pc), corresponding to 
a parallax smaller than approximately 5.5 mas. For that reason, we 
can exclude three of the YLA - Tuc-Hor, /3 Pic-Cap and AB Dor - 
for which the current distance to RBS 1223 would be smaller than 
100 pc. 8 From the r-dmin contour plot for HD 141569, a potential 
supernova in HD 141569 would have occurred near the edge of the 
group our even outside (between 5 and 21 pc, the radius of HD 
141569 is Ki 16 pc) and thus is unlikely. 



8 As for RX J1605.3+3249, the distance of RBS 1223 is too poorly known 
to estimate probabilities for the potential parent associations. For that rea- 
son, we give an individual discussion. 
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Table 10. Potential parent associations of RBS 1223, columns as in Table 4. 



Association d min r v r ju* fj,s v apaC e t cIq a 5 

[pc] [Myr] [km/s] [mas/yr] [mas/yr] [km/s] [mas] [pc] [°] [°] 

Tuc-Hor 23... 30 0.17... 0.31 463+^ -206 ± 19 84 ± 18 473±gg° 10.85±|"_f| 37l 4 339.68lg 4 6 g 2 -29. 201^3 

P Pic-Cap 5... 10 0.13... 0.22 482lg° 6 -206 ± 19 82 ± 19 490lgg° i2.O6l2.25 19 -3 S 16 - 23 !^ 1 -29.63l3;g°, 

HD 141569 5... 21 0.27... 0.40 332lg° 3 -209 ± 19 85 ± 20 393lgg° 5.09l°,;g 8 84... 122 236.46l g '; 7 )g ) 1.55±Hl 

Ext. RCrA 16... 40 0.38... 0.58 345lgg 4 -209 ± 20 85 ± 18 39ll 103 " 5.78±\\H 80l^ 294.821 7 ,; 24 -22.941 4 ; 33 

ABDor 15... 28 0.15... 0.24 437±g5 -206 ± 20 83 ± 19 446lg° 4 ll.56l2.S5 18 -i 225 . . . 358 -28.86lg 1 8 g 7 

US 42... 64 0.38... 0.55 35lli°4 -210 ± 20 83 ± 18 420l^s 4m ~i.27 1091^ 254.09I2.34 - 6 - 29 -tol 



Table 11. Results for RBS 1223 and Set OB2A and B. Minimum separations have been constrained to d m in $J 50 pc. Column 2 gives the smallest separation 
dmin found after 2 million runs, column 3 gives the maximum of the time (r) histogram along with its 68% interval (note that this is not a la error). Columns 4 
to 10 indicate the current neutron star properties needed to yield small separations to the particular association as well as the position of the potential supernova 
(see also Table 4). Note that receding radial velocities are positive contrary to Motch et al. (2009). 
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Indeed, the case Ext. R CrA is very promising, since the shape 
of the d min distribution is in good agreement with a three- 
dimensional Gaussian distribution with /1 = and a — 13.5 (see 
Fig. 10(a), equation 2). Hence, a potential supernova which created 
RBS 1223 would have occurred near the centre of the association. 
Assuming this, RBS 1223 would now be « 0.5 Myr old. Its 
current radial velocity would be w 350 km/s (space velocity 
« 400 km/s). Furthermore, it is very likely that Ext. R CrA 
experienced a supernova in the near past owing to the existence of 
B-type stars (highest content within the YLA, F08). 
In Table 10 (bottom) we also give the results for the US association 
which has also been suggested as a possible birth place for RBS 
1223 by Motch et al. (2009). Although the d min histogram in 
this case does not support a birth in US (a supernova is predicted 
at a distance from the US centre of 51 ± 12 pc, see Fig. 10(b)), 
it seems possible that the neutron star was born in the Sco-Cen 
complex (similar to RX J1605. 3+3249, see previous section) in 
the vicinity of US. In this case, the present distance of RBS 1223 
to the Sun would be slightly larger as for Ext. R CrA (w 215 pc) 
whereas age and radial velocity are similar to those for Ext. R CrA 
(r « 0.5 Myr, v r w 350 km/s, space velocity « 420 km/s). 
Adopting a near birth scenario either in the Sco-Cen complex or 
Ext. R CrA 0.5 Myr ago (Fig. 1 1), we can estimate the mass of the 
progenitor of RBS 1223. Given association ages of 5 to 10 Myr for 
US and 10 to 15 Myr for Ext. R CrA, the mass would be 17 (B0V) 
to 45 M Q (05.5V) (US) or 12 (B1V) to 30 M Q (07V) (Ext. R 
CrA) (T80, MM89, K97; spectral types from Schmidt-Kaler 1982). 

However, since the distance to RBS 1223 is very uncertain 
and following P07 can only hardly get constrained through uh, 
we eventually cannot completely exclude the associations in the 
lower part of Table 9. Infact, in the case RBS 1223 P07 reached the 
limit of their modelling which was set at 1 kpc. Thus, the neutron 
star might be located at somewhat larger distances which would be 
consistent with its origin lying in the Set region. If so, we can again 
estimate the mass of the progenitor star adopting an age of RBS 
1223 of 1 Myr (Fig. 1 1) and Set OB2 of 6 Myr (Schild & Maeder 
1985). Models from T80, MM89 and K97 yield masses of 44 M , 
38 M and 32 M Q (^06 on the main sequence). 



5.5 Remarks 

It has to be mentioned that any association identified for which it 
is possible to trace our neutron stars back such that they reach the 
association boundary (although it might be only fulfilled for a very 
small fraction of runs) is a potential parent association. However, 
for a few associations, theoretical curves (equations 1 and 2) fit the 
d m in histograms remarkably well. Hence, we conclude that they 
reflect the real separation at the encounter time (flight time). 
Furthermore, it is important to bear in mind that the initially input 
distribution for radial velocities (Ho05) somewhat biases the out- 
coming radial velocities to median ranged values (w ±400 km/s, 
depending on the transverse velocity). Due to this bias, the range of 
t (flight time) might be underestimated. Nonetheless, our method 
allows us to also identify associations requiring either very small 
or high radial velocities. 

Compared to previous investigations of the origin of M7 
members as done by e.g. Walter & Lattimer (2002), Motch 
et al. (2003), Kaplan et al. (2007) and Motch et al. (2009), our 
calculations do not only cover certain radial velocity ranges, but 
the whole spectrum of possible radial velocities ranging from 
— 1, 500 to +1, 500 km/s. To account for the velocity distribution 
for neutron stars we used the probability distribution of Ho05. 
Furthermore, we vary the distances for RX J1605.3+3249 and RBS 
1223 in a wide range to account for their large uncertainties. 
We compiled a list of 140 OB associations and clusters to cover 
most of the potential birth sites of neutron stars. 

For estimating the supernova progenitor star mass from the 
difference between association age and kinematic age, we assumed 
that all stars formed about at the same time in that association. This 
assumption is justfied, because once an O or B star is formed, its 
strong wind and slightly later its supernova shock front will blow 
away all the remaining gas, so that star formation will cease in that 
region. This assumption holds for OB associations and is consis- 



9 However, some neutron stars may have been former massive runaway 
stars. 
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tent with observations. It does not neccessary hold for associations 
without an O- or B-type star, like possibly TWA, which however 
may have had an early B-type star. 



6 CONCLUSIONS 

We reviewed the prominent neutron star/ runaway case PSR 
B 1929+10/ £ Oph (HOI) and reinvestigated it with more recent 
data for the pulsar. Although we did not put constraints on the 
radial velocity, we could confirm previous results from Bobylev 
(2008) and eventually derive a radial velocity of w 250km/s 
of the pulsar which is needed to find the two objects at nearly 
the same position in US at the same time in the past. However, 
we should bear in mind that we had to increase the errors of the 
parallax and proper motion components. If they are truly as small 
as published by Chatterjee et al. (2004) the scenario is less likely. 

We investigated the motion of four members of the M7 
and for each neutron star simultaneously the motion of 140 OB 
associations and clusters to confine potential parent associations of 
the neutron stars. A summary of our results is given in Table 12 
(for individual discussions please see previous sections). 
We find that RX J 1856.5-3754 most probably originated from US 
about 0.3 Myr ago. RX J0720.4-3 125 was very likely born in TWA 
about 0.4 Myr ago or may come from Tr 10 where the supernova 
then would have occured 0.5 Myr in the past. RX J1605. 3+3249 
and RBS 1223 also were seemingly born from a close young 
association such as the Scorpius-Centraurus complex or Ext. R 
CrA. For RBS 1223 also a birth in Scutum OB2 is possible. 

Thus, we find that the most probable birth places were lo- 
cated close to the Sun. The supernova distances are consistent with 
the Local Bubble and hence may have contributed to its forma- 
tion or reheating (Berghofer & Breitschwerdt 2002; F08). There 
is no doubt that the Sco OB2 association to which US belongs 
experienced some supernovae in the past (Mafz-Apellaniz 2001; 
Berghofer & Breitschwerdt 2002; F08). There have been other in- 
vestigations showing that many neutron stars may have been born 
within the Gould Belt (Popov et al. 2003), a torus-like structure 
around the Sun with a radius of w 600 pc (Torra et al. 2000), to 
which most of the associations just mentioned also belong. 
Some of the four investigated M7 neutron stars may have formed 
within « 250 pc within the last 0.5 Myr, so that they could have 
contributed to the 1() Be and 60 Fe found in the Earth's crust. 

Adopting a birth scenario of each neutron star, we can 
derive estimations of the mass of the progenitor star assuming 
contemporary star formation in the association as well as the 
kinematic neutron star age which we here always find to be lower 
than the characteristic spin-down age (Table 13) that is an upper 
limit to the true age. From cooling models we know that, given 
the surface temperature, for the neutron stars investigated here, 
the characteristic age is too high. This is also plausible since those 
objects are still very young. Thus, our kinematic ages better fit 
with cooling models (see Fig. 12). 

Since RX 11856.5-3754 is the coolest of the four M7 members 
investigated (RX J0720.4-3125, RX J1605.3+3249 and RBS 1223 
show similar temperatures), it should be the oldest. However, it 
seems, that RX J1856.5-3754 is slightly younger than the other 
three (or has a similar kinematic age if we adopt the parallax of 
Walter et al. 2000). This might be owing to different parameters at 
birth or different cooling. Also the kinematic age of RX J1856.5- 



Table 13. Kinematic ages r kin of the four M7 members compared with 
their characteristic spin-down ages r char (taken from the ATNF pulsar 
database 10 , Manchester et al. 2005). The fifth column gives the effective 
temperature as listed in Table 1 of Haberl (2005) (see references therein). 
Column 6 gives the median masses (for TWA a mass consistent with its 
mass function) of the progenitor stars predicted using evolutionary models 
from T80, MM89 and K97 (see text) given the difference between the ages 
of the associations and the time since the predicted supernovae. 
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3754 is of about a factor of 13 smaller than its characteristic age 
which is a large difference. This may indicate further doubts on 
the characteristic age. Mori & Ruderman (2003) propose that 
RX J1856. 5-3754 is an isolated magnetar that has spun down by 
the propeller effect. In this case, dipole braking, which is assumed 
to calculate characteristic ages, is no longer valid. Moreover, 
RX J1856. 5-3754 does only show a very small pulse fraction 
making it difficult to derive P (van Kerkwijk & Kaplan 2008). 
In addition, effective temperatures may be influenced by hot spots 
and do not reflect the real surface temperature. Infact, this is 
the case for RX J0720.4-3125 (Haberl et al. 2006; Hohle et al. 
2009) and RBS 1223 (Schwope et al. 2005; Haberl 2007) which 
appear hotter, thus must be actually shifted somewhat to smaller 
temperatures (even more consistent with cooling curves). 
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Table 12. Summary of the results discussed in section 5. Given here are those associations that we conclude to be the most probable birth places of the four M7 
members. Columns 3 to 5 give the position of the potential supernova (SN) as equatorial coordinates and distance to the Sun, column 6 gives the approximate 
time in the past at which the SN may have occurred and columns 7 to 10 give constraints on the neutron star parameters as they would be needed to reach the 
respective association at the given time. 
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Figure 12. The four M7 members inserted into a cooling diagram. Filled 
stars mark the characteristic spin-down age (see Table 13) whereas hor- 
izontal lines characterise an area of the kinematic age (lower and upper 
values from associations in tables of section 5). Open diamonds show the 
kinematic age for the associations summarised in Table 12. Effective tem- 
peratures can be found in Table 13. 

The purple set of cooling curves was adopted from Popov et al. (2006) (solid 
lines, for masses of 1.05, 1.13, 1.22, 1.28, 1.35, 1.45, 1.55, 1.65 and 
1.75 Mq from top to bottom; Model from Grigorian et al. 2005), the green 
set has been kindly provided by A. D. Kaminker (dashed lines, includes su- 
perconductive protons and neutrons, for masses from 1.1 to 1.9 Mq from 
top to bottom; see also Gusakov et al. 2005). 

(Manchester et al. 2005). 
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Table Al. The sample of OB associations and clusters. The table gives the consecutive number, the designation, the position on the sky (J2000 galactic 
longitude I and latitude ft), the parallax n, heliocentric velocity components U, V and W as well as estimated ages and spatial extensions either as found in 
the literature or obtained from angular extensions and distances. Besides Sco OB2 and YLA, U, V and W were computed as in Johnson & Soderblom (1987) 
(errors include proper motion and radial velocity uncertainties only). Note that ages are often uncertain to a factor up to two. Many OB associations contain 
several clusters of different ages (many of them listed here as well). 

The last column lists the references used for each group: A06 - Aharonian et al. (2006), B08 - Bobylev (2008), BaL95 - Balona & Laney (1995), BH89 - 
Blaha & Humphreys (1989) (according to Dambis et al. (2001) distances taken from this publication have been reduced by 20%), Bh07 - Bhavya et al. (2007), 
B156 - Blaauw (1956), B191 - Blaauw (1991), BL04 - Beshenov & Loktin (2004), Br99 - Brown et al. (1999), BYN06 - Barrado Y Navascues (2006), C06 - 
Carraro et al. (2006), CH08 - Carpenter & Hodapp (2008), C174 - Claria (1974), Co02 - Contreras et al. (2002), D01 - Dambis et al. (2001), dG89 - de Geus 
et al. (1989), Di02 - Dias et al. (2002), dlR06 - de la Reza et al. (2006), dZ99 - de Zeeuw et al. (1999), E07 - Eisenbeiss (2007), F08 - Fernandez et al. (2008), 
Fu04 - Fuhrmann (2004), G94 - Garcia (1994), GB08 - Gvaramadze & Bomans (2008), GH08 - Gregorio-Hetem (2008), Gr71 - Graham (1971), GS92 - 
Garmany & Stencel (1992), H01 - Hoogerwerf et al. (2001), Ha72 - Havlen (1972), HR76 - Herbst & Racine (1976), Hu78 - Humphreys (1978), J05 - Jilinski 
et al. (2005), K05 - Kharchenko et al. (2005), K07 - Kharchenko et al. (2007), KG94 - Kaltcheva & Georgiev (1994), KM07 - Kaltcheva & Makarov (2007), 
L06 - Lopez-Santiago et al. (2006), LB01 - Loktin & Beshenov (2001), LB03 - Loktin & Beshenov (2003), Lo06 - Lodieu et al. (2006), Loz86 - Lozinskaia 
et al. (1986), LS04 - Luhman & Steeghs (2004), Lu08 - Luhman (2008), M02 - Madsen et al. (2002), Ma07 - Makarov (2007), Mae87 - Maeder (1987), 
Mas95 - Massey et al. (1995), MasOl - Massey et al. (2001), MC68 - MacConnell (1968), ME95 - Mel'nik & Efremov (1995), Me07 - Menten et al. (2007), 
Mer08 - Mermilliod et al. (2008), Meri04 - Merin et al. (2004), MP03 - Mermilliod & Paunzen (2003, WEBDA database, http://www.univie.ac.at/webda), 
MV73 - Moffat & Vogt (1973), Mo53 - Morgan et al. (1953), N00 - Neuhauser et al. (2000), P91 - Perry et al. (1991), Pr02 - Preibisch et al. (2002), S06 

- Smith (2006), Sa03 - Sartori et al. (2003), SM85 - Schild & Maeder (1985), SN00 - Stelzer & Neuhauser (2000), So04 - Southworth et al. (2004), Sod98 

- Soderblom et al. (1998), T76 - Turner (1976), T79 - Turner (1979), T80a - Turner (1980), T80b - Turner et al. (1980), E07 - Eisenbeiss (2007), Te99 - 
Terranegra et al. (1999), U01 - Uyamker et al. (2001), W07 - Wolff et al. (2007), We63 - Westerlund (1963), Web99 - Webb et al. (1999), Wh97 - Whittet 
et al. (1997), Wo08 - Wolk et al. (2008), ZOla - Zuckerman et al. (2001), ZOlb - Zuckerman et al. (2001). 



Nr 


Name 


![°] 


M°] 


7r [mas] 


U [km 


/*] 


V [km/s] 


W [km/s] 


Age [Myr] 


®[pc] 


Ref. 


Scorpius OB2 


1 


US 


351.07 


19.43 


6.65 


-6.7± 


5.9 


-16.0 ± 3.5 


-8.0± 2.7 


5-10 


30 


dG89, dZ99, Br99, 
























M02, Pr02, Sa03, F08 


2 


UCL 


330.51 


12.86 


6.96 


-6.8± 


4.6 


-19.3 ± 4.7 


-5.7 ± 2.5 


10-20 


65 


dG89, dZ99, Br99, 
























M02, Sa03, F08 


3 


LCC 


301.54 


6.74 


8.38 


-8.2± 


5.1 


-18.6 ± 7.3 


-6.4± 2.6 


13-20 


45 


dG89, dZ99, Br99, 
























M02, Sa03, F08 


Young Local Associations 


4 


TWA 


291.61 


20.22 


16.46 


-9.7± 


4.1 


-17.1 ± 3.1 


-4.8± 3.7 


8-20 


66 


Sod98, Web99, dlR06, 
























BYN06, F08 


5 


Tuc-Hor 


296.57 


-51.72 


23.09 


-10.1 ± 


2.4 


-20.7 ± 2.3 


-2.5 ± 3.8 


10-40 


100 


SN00, ZOlb, Ma07, 
























F08 


6 


/3 Pic-Cap 


330.95 


-55.54 


54.97 


-10. 8± 


3.4 


-15.9 ± 1.2 


-9.8 ± 2.5 


8-34 


113 


ZOla, Ma07, F08 


7 


e Cha 


300.43 


-15.08 


10.41 


-8.6± 


3.6 


-18.6 ± 0.8 


-9.3± 1.7 


5-15 


55 


Te99, J05, F08 


8 


rj Cha 


292.42 


-21.45 


10.76 


-12.2 ± 


0.0 


-18.1 ± 0.9 


-10. 1± 0.5 


5-8 


13 


LS04, J05, F08 


9 


HD 141569 


7.40 


39.50 


9.94 


-5.4± 


1.5 


-15.6 ± 2.6 


-4.4± 0.8 


5 


31 


Meri04, F08 


10 


Ext. R CrA 


359.41 


-17.19 


9.85 


-0.1± 


6.4 


-14.8 ± 1.4 


-10. 1± 3.3 


10-15 


62 


N00, F08 


11 


AB Dor 


146.31 


-59.00 


71.43 


-7.4± 


3.2 


-27. 4± 3.2 


-12. 9± 6.4 


30-120 


85 


L06, Ma07, F08 


Hercules-Lyrae association 


12 


Her-Lyr 


180.00 


68.20 


52.10 


-13. 6± 


3.7 


-22. 7± 3.4 


-7.9± 5.2 


40-200 


26 


Fu04, L06, E07 


Other OB associations and clusters 


13 


Sgr OB5 


0.00 


-1.19 


0.41 


-6.5 ±10.8 


1.2 ±11.6 


1.3±15.0 


6-12 


221 


SM85, BH89, D01 


14 


NGC 6530 


6.07 


-1.33 


0.76 


-12.1 ± 


6.3 


-16.2 ± 1.0 


-14. 3± 0.8 


4-5 


11 


LB03, K05 


15 


Sgr OBI 


7.58 


-0.78 


0.65 


-8.0 ± 


1.8 


-12. 8± 1.5 


-9.5 ± 1.5 


5-8 


207 


SM85, BH89, D01 


16 


Sgr OB7 


10.71 


-1.52 


0.54 


-7.4± 


3.1 


-3.2 ± 2.3 


-18. 6± 2.4 


4-5 


58 


K05 


17 


Sgr OB4 


12.10 


-0.99 


0.52 


1.7 ± 


3.6 


-6.2 ± 7.2 


-7.3 ±12.8 


<10 j 


54 


Hu78, BH89, D01 


18 


Sgr OB6 a 


14.18 


1.28 


0.62 


-10.8 ± 


3.9 


-2.8 ± 7.4 


-16.0 ± 8.6 


<10-> 


34 


Hu78, BH89, ME95 


19 


M17 


15.04 


-0.68 


0.55 


-20.5 ±12.7 


9.9 ± 7.4 


-21. 6± 8.9 


6 


19 


K05 


20 


Set OBI 


16.73 


0.00 


1.88 


-4.3± 


4.9 


-3.1 ± 1.6 


-2.1± 1.3 


8-13J 


71 


Hu78, BH89, D01 


21 


NGC 6611 


16.95 


0.79 


0.58 


17. 2± 


0.6 


5.5 ± 0.7 


-5.6± 0.8 


1-5 


13 


Mas95 LB03, K05, 
























GB08 


22 


Set OB3 


17.32 


-0.84 


0.74 


7.4 ± 


6.1 


-3.7± 2.6 


-3.9± 3.8 


4.5 j 


33 


Hu78, BH89, D01 


23 


Ser OB2 


18.23 


1.66 


0.63 


-1.6± 


1.8 


-6.9 ± 2.2 


-10.7 ± 0.8 


4.5-i 


61 


Hu78, BH89, D01 


24 


NGC 6604 


18.25 


1.71 


0.59 


11.9 ± 


7.1 


-5.2± 2.7 


1.5 ± 1.5 


4-5 


15 


LB03, K05 


25 


Set OB2 


23.21 


-0.53 


0.63 


-8.7± 


7.5 


-7.8± 8.1 


-4.4± 3.8 


<6 


61 


SM85, BH89, D01 


26 


Tr35 


28.28 


0.01 


0.83 


4.1 ± 


1.9 


-17.6 ± 3.3 


-1.0 ± 3.6 


73 


2 


Di02 


27 


Col 359 


29.66 


12.72 


2.22 


6.5± 


0.3 


-15.8 ± 0.5 


-10.8 ± 0.7 


28-80 


17 


LB03, K05, K07, 
























Lo06, B08 


28 


IC 4665 


30.62 


17.11 


2.60 


2.0± 


1.5 


-18.7 ± 0.9 


-9.7 ± 0.6 


43 


13 


H01.LB03, K05 


29 


Vul OBI 


59.41 


-0.12 


0.63 


34.4 ± 


4.8 


-11.3 ± 5.4 


-0.8 ± 1.5 


10-16 j 


230 


Hu78, T80a, BH89, 
























D01 


30 


NGC 6823 


59.40 


-0.15 


0.43 


63.9 ± 


3.0 


-18.0 ± 1.8 


-0.8± 4.1 


2-7 


68 


K05, W07 


31 


Vul OB4 


60.66 


-0.17 


1.25 


12.1 ± 


2.9 


-10.2 ± 3.9 


-7.2 ± 2.7 


10 


35 


T80a, BH89, D01 


32 


Cyg OB3 


72.81 


2.00 


0.43 


77.6 ± 


2.2 


-32. 9± 2.1 


-13. 6± 3.3 


8-12 


106 


BH89,D01,U01,So04 
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Table Al: - Continued. - 



Nr 


Name 


in 


b[°] 


7r [mas] 


U [km/s] 


V [km/s] 


W [km/s] 


Age [Myr] 


<Z>[pc] 


Ref. 


33 


NGC 6871 


72.65 


2.05 


0.64 


49. 5± 1.9 


— 25. 8± 2.2 


— 10.4 ± 2.0 


2-10 


13 


Mas95, LB03, So04, 






















K05 


34 


Byurakan 1 


72.75 


1.75 


0.63 


48. 9± 3.8 


-24.0 zt 3.0 


-3.4 ± 3.4 


18 


12 


K05 


35 


Byurakan 2 


72.75 


1.34 


0.58 


49.8 ± 4.3 


-38.4 ± 9.1 


-3.3 ± 3.2 


5 


11 


K05, Bh07 


36 


NGC 6883 


73.28 


1.18 


0.48 


68.8 ± 3.0 


-31.3 =t 4.2 


-4.6 ± 2.8 


34 


8 


K05, A06 


37 


Cyg OBI 


75.90 


1.12 


0.59 


45.4 ± 2.4 
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Table A 1 : - Continued. - 



Nr Name l[°] b[°] n [mas] U [km/s] V [km/s] W [km/s] Age [Myr] ®[pc] Ref. 
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.9 


2 


6± 1 


9 


-15.3 ± 


1 
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a v r from ME95, proper motion components derived from Hu78 members as median ± average deviation from median. 



Proper motion components and v r derived from W63 members as median ± average deviation from median. 
c Coordinates and parallax denote the mean of Cas-Tau members given by dZ99, velocity components reflect their median U,V,W. 
d Proper motion components and v r derived from CH08 members as median ± average deviation from median. 
G Velocity components were derived from v r only. 

f Proper motion components and v r derived from Hu78 members as median ± average deviation from median. 
s Proper motion components and v r derived from Wh97 members as median ± average deviation from median. 
h v T derived from Webda (MP03) members as median ± average deviation from median. 

1 Proper motion components and v r derived from MasOl members; proper motion (2 stars): mean ± maximum deviation, v T : for only one star available (without deviation). 

J We derived ages from HR diagrams from memberlists either from Hu78 or GS92 (Mo53 for Cyg OB4, Gr7 1 for Mon OB3) by fitting a set of theoretical Padova isochrones (Girardi et al., 2002) assuming solar 

metallicity. The fitting algorithm is based on the least-squares method and uses stellar magnitudes as weights. 

k Ori OB la: 11 Myr, Ori OB lb: 2 Myr, OriOBlc: 5 — 11 Myr, Ori OBld (= Trapezium cluster): S3 1 Myr. 



APPENDIX A: SAMPLE OF OB ASSOCIATIONS AND 
CLUSTERS 

For our investigation we selected OB associations and clusters 
within 3kpc from the Sun (Table Al; for detailes on the selection 
criteria, please see section 2). 



APPENDIX B: DERIVATION OF PRESENT-DAY 
NEUTRON STAR PARAMETERS AND SUPERNOVA 
POSITION - EXAMPLE RX J1856.5-3754 COMING FROM 
US 

In Tables 4, 6, 8 and 10 we list the current proper motions (/x* , 
Us ; columns 5 and 6), radial velocities (v r ; column 4) and paral- 
laxes (7r; column 8) which the neutron stars would have if they 
originated from the particular associations as well as the positions 
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Figure Bl. Distribution of the number of runs in the T-d m i n space (ex- 
cluding runs yielding r = Q) for RX J1856. 5-3754 and US. The ellipse 
marks the region for determination of the parameters given in Table 4 and 
Fig. B2. 



of the predicted supernova (distance to the Sun and equatorial 
coordinates; columns 9 to 11). Here we describe how those values 
are obtained. 

We first define an area within the T-d m in contour plot such that 
its boundaries reflect a 68% decline from the peak (columns 2 and 
3 of the respective tables). Parameter values in columns 4 to 11 
were derived by selecting the input parameters (parallax, proper 
motion, radial velocity) which correspond to this defined region. 
The distance to the Sun as well as equatorial coordinates at the 
time of the supernova were then calculated for each parameter set. 
From the histogram of each parameter we obtain its value and error 
by drawing an interpolation curve to better characterise the shape 
of the distribution. The "mean" of the parameter is then given by 
the maximum of the curve. The error intervals include about 68% 
of the histogram area (note that these are not la errors). In the case 
of multiple peaks in the histogram, the interval between the two 
values which reflect a 68% decrease of the maximum is given. 
Fig. B 1 gives an example of a T-d m i n contour diagram indicating 
the area described above. The plot shows the distribution of mini- 
mum separations between RX Jl 856.5-3754 and the US centre in 
the T-d m in space. Histograms corresponding to the entries for US 
in columns 4 to 6 and 8 to 1 1 of Table 4 are shown in Fig. B2. 

The probability of an association being the birth place of a 
neutron star is roughly estimated from the probability of the paral- 
lax value needed for the scenario (Gaussian distribution according 
to the literature value) and the probability of the necessary space 
velocity (adopting the distribution by Ho05). Furthermore, we are 
taking into account the distance of the potential supernova from the 
centre of the association since due to a decrease of the star density 
with increasing distance from the centre of an association, super- 
novae occur more likely near the centre. To account for that we 
make a rough estimate of the probability of a potential supernova 
position by taking a Gaussian distribution with a being half the ra- 
dius of the association (weighted according to the number of mem- 
ber stars). This is justified since radii of associations are derived 
from star density profiles in most cases. The three probabilities are 
then multiplied. 
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Figure B2. (a)-(d): Distributions of present-day parameters for RX J1856.5-3754 supporting that it was created in a supernova in US about 0.3 Myr ago. 
(e)-(g): Position of the potential supernova. Lines are interpolations to easier see the shape of the histogram and determine the confidence intervals of the 
parameters. 



